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We consider a tunnel junction formed between a fixed electrode and an oscillating one. Accu- 
mulation of the charge on the junction capacitor induces a force on the nano-mechanical oscillator. 
The junction is voltage biased and connected in series with an impedance Z(u). We discuss how the 
picture of Coulomb blockade is modified by the presence of the oscillator. Quantum fluctuations of 
the mechanical oscillator modify the I-V characteristics particularly in the strong Coulomb blockade 
limit. We show that the oscillator can be taken into account by a simple modification of the effective 
impedance of the circuit. We discuss in some details the case of a single inductance Z(lj) = iLu and 
of a constant resistance Z(u>) = R. With little modifications the theory applies also to incoherent 
transport in Josephson junctions in the tunneling limit. 

PACS numbers: 73.23.-b, 85.85.+J, 74.50.+r 



I. INTRODUCTION 



The improved ability to build electronic devices on 
the nanometer scale opens the perspective to exploit and 
study the coupling between electronic transport and me- 
chanical degrees of freedom!^ The elementary system 
typically considered in nano-clcctromcchanics is the har- 
monic oscillator coupled in some way to the electronic 
degrees of freedom. A widely studied device is the single- 
electron transistor which consists in a small metallic is- 
land connected to the leads through two tunnel barriers. 
The oscillating part of the device can be cither a nearby 
gate electrode, or the island itself. The motion of the 
device modifies mainly two quantities: the capacitances 
and the (bare) tunneling rates. According to which de- 
pendence dominates, several effects have been predicted 
or observed. When only the capacitance is modified by 
the oscillation, it has been shown in the weak coupling 
limit that the variance of the oscillator position satis- 
fies the equipartition theorem with the voltage bias re- 
placing the temperature^ Under specific conditions, the 
system may also undergo dynamical instabilities! 4 -!^ In 
the strong coupling limit, a new kind of current block- 
ade (Frank-Condon blockade) has been predictedf 7 -^^ 
This effect also persists in the classical limit] 11 ' 12 ' 13 When 
the oscillation modifies the distance between the metal- 
lic leads, it has been shown that a dynamical instability 



can occur 
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This instability (called shuttle 
instability) is characterized by a synchronization of the 
electronic transport with the mechanical oscillations, as 
is clearly seen from the full counting statistics of charge 
transport . 21 ' 22 Observation of shuttling is difficult, since 
the distance between the metallic leads must be com- 
parable with the tunneling length in order to produce 



a measurable current] 23 ' 24 ' 25 ' 26 Very recently indications 
of shuttling have been found by the authors of Ref. 
probably in an intermediate regime between tunneling 
and field emission. The wide interest for the single elec- 
tron transistor comes from its high sensitivity to a small 
variation of the gate voltage. Thus, it can be used, for in- 
stance, as a sensitive displacement detector— and it has 
been shown to reach the quantum back action limits 

An even simpler device is the single tunnel junc- 
tion where the distance between the two metals con- 
stitutes the mechanical degree of freedom. The effect 
of the position dependence of the tunneling resistance 
has been already considered in the literature ; 30 ' 31 ' 32 both 
for the current and the noise of the device. In this 
case, the force acting on the oscillator originates from 
the electron-momentum transfer. Recently, the position 
fluctuation has been detected by measuring the current 
fluctuations.— This experiment probed the back action of 
the current crossing the tunnel junction on the oscillator. 
Apparently, the position dependence of the tunneling ma- 
trix elements is not sufficient to explain the intensity of 
the back action. To our knowledge, the effect of the force 
coming from the variation of the capacitance C as a func- 
tion of the distance between the leads has not been con- 
sidered in this context so far. A particularly interesting 
case is when the device is voltage biased in series with 
an impedance Z(uS) leading, under certain conditions, to 
the Coulomb blockade physics^ 4 - Then, similarly to what 
happens for the single electron transistor, one can expect 
that a small position dependence of the capacitance can 
affect the current- voltage characteristics of the device. 

In this paper, we consider the effect of the mechani- 
cal oscillator on the I-V curve. We will work within the 
assumptions of the standard Coulomb blockade theory^ 



2 



Specifically, we assume that the tunneling resistance is 
sufficiently large so that the oscillator and the electro- 
magnetic environment have time to relax to equilibrium 
between two tunneling events. This is the opposite limit 
with respect to the one considered in Refs. I3lll32l for the 
tunneling dependent Hamiltonian. There, the tunneling 
is so frequent that the system may reach a stationary 
off-equilibrium state. Restricting to the linear coupling 
of the oscillator with the electromagnetic modes, we will 
discuss the effect of the quantum fluctuations of the me- 
chanical oscillator on the current. 

The paper is organized as follows. In Sect. [IT] we de- 
scribe the model and its range of validity. We show that 
within a linear approximation the mechanical oscillator 
can be accounted by an effective impedance. In Sect. IIIII 
we consider two specific cases of Z(io) and discuss the 
resulting conductance at zero temperature. In Sect. IIVI 
we consider the case of an oscillating Josephson junction 
by extending the results of previous sections to the su- 
perconducting case. Section fVl gives our conclusions. 



II. MODEL 

We consider the electronic transport through a circuit 
composed of a tunnel junction in series with an arbitrary 
impedance Z(ui) and biased with an ideal voltage source 
V. We assume that one of the two metallic leads forming 
the tunnel junction is free to oscillate and can thus mod- 
ify its distance with respect the other one. The single 
degree of freedom describing the oscillation is given by 
x, the oscillator position. Its effective mass and spring 
constant are m and k, respectively. The tunneling sur- 
face is very small. Thus, it is important to account for 
the effective capacitor forming at the tunnel barrier. Its 
capacitance C{x) depends on the position of the oscilla- 
tor. The tunneling resistance Rt{x) may also depend on 
x. In the following we will neglect this dependence by 
showing that its effect is negligible with respect to that 
of C(x) in the regime of rare tunneling events that we 
are considering. The device is depicted schematically in 
Fig.fTJ 



A. Time scales 

As anticipated in the Introduction we will work within 
the same hypothesis that are commonly used to de- 
scribe Coulomb blockade. Specifically, we assume that 
the relaxation times for the electromagnetic environment, 
^-/iRC = RC, an d for the oscillator, I/7, are shorter 
than the average time between two tunneling events, 
tt ~ eRr/V. Here, R ~ Re[Z] is the typical resistance 
associated with the impedance Z and 7 is the damp- 
ing coe fficient of the mechanical oscillator with frequency 
uj m = sjk/m. If I/jrc, I/7 <C Vp, the system has time 
to come back to the equilibrium state after an electron 
has crossed the tunnel barrier. This allows us to calculate 




FIG. 1: Schematic representation of the system. A tunnel- 
ing contact is formed between an oscillating electrode (here 
depicted as a STM tip) and a fixed surface. The oscillator 
is modelled by a spring constant k and an effective mass m. 
The junction between the tip and the plate is characterized 
by an effective capacitance C(x) and resistance Rt(x), where 
x denotes the position of the tip. The circuit is closed by an 
impedance Z and a voltage source V . 

the average current on the equilibrium quantum state of 
the environment. It gives the following limitations on 
the voltage: eV <C Ec(Rt/R) and eV <C e 2 7-fi/r, where 
Eq = e 2 /2C is the Coulomb energy. Since we are inter- 
ested in the range of voltages around Eq /e we assume 
that Rt ^> R and Rt S> l/Cty The first condition 
is the standard tunneling condition, the second involves 
the mechanical damping. One should remember that the 
damping of the oscillator also depends on the impedance 
of the circuit. Thus, even if the oscillator is not damped 
by a mechanical source, its quality factor would not be 
infinite. In the limit of large Rt the above conditions 
can thus always be satisfied. 

Let us now discuss the effect of the position dependence 
of the tunneling resistance. To evaluate the importance 
of this effect for the dynamics of the junction we consider 
the effective friction coefficient generated by the position 
dependence of Rt'M- 

where Rq = 2irh/ e 2 is the quantum of resistance, x 2 = 
h/2muj m is the (square of the) zero point motion am- 
plitude of the oscillator and At is the tunneling length 
defined by A^ 1 = {cLRt / dx) / Rt ■ The explicit expression 
of the friction coefficient induced by the electromagnetic 
dissipation is given below for a specific case [cf. Eq. (J25}]. 
Independently on its explicit form, what matters is that 
it does not depend on the tunneling resistance since the 
damping is due to the environment and the capacitivc 
coupling. Thus the ratio 77 1 /"f can always be made small 
by choosing a large enough tunneling resistance. Then 
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we may safely neglect the position dependence of the tun- 
neling matrix elements. We will come back to this point 
at the end of Section ITTT1 



B. Derivation of the general expression for the 
current 

The Hamiltonian describing the system is the follow- 
ing: 



H = Hi 



H n 



H 



EM- 



(2) 



where H 
A 



,i<t>(t)J c ^ 



Ht with Ht = 
Here, k and q label 



3 ~ ^k ' 

At and A = J2 k q t kq e* 
the femionic eigenstates in the fixed and oscillating elec- 
trode, respectively, c k is a destruction operator in state 
k, with energy and t kq the tunneling matrix clement. 
The phase <jj(t) is obtained by integrating the voltage dif- 
ference at the junction: <fi(t) — J_ (eVj(t')/h)dt'. For 
a constant voltage biased junction <j>(t) = eVt/h. How- 
ever, due to the impedance in series with the junction, 
the voltage fluctuates and, in general, Vj is different from 
V. The Coulomb blockade physics stems from the behav- 
ior of the phase <fi and, specifically, from its correlation 
function. The current through the junction can be ob- 
tained by perturbation theory in Ht- The operator for 
the current from the mobile electrode to the fixed one 
reads / = — ie(A — A')/h. Linear response theory gives: 



7 = 4 Em 2 / dt 

11 kg J -°° 



-*<*>e*<°>) /*(! - /,) - (eW^W) f q (l - f k ) 



(3) 



where f k . q = l/(e Kk " + 1) are the Fermi distributions in 
the leads, and (3 is the inverse temperature. In order to 
include the quantum fluctuation of the electromagnetic 
fields, the phases in this expression must be regarded as 
operators whose dynamics is determined by the rest of 
the Hamiltonian We will restrict to the quadratic 
part of the Hamiltonian. Within this approximation, the 
fields 4> satisfy Wick theorem and expression can be 
simplified by making use of the exact relation: 



i<£(0)\ 



,J(t)-J(0) 



(4) 



where (j) = <j>{t) - eVt/H and J(t) = U(t)^)(0) 

The current can then be expressed in terms of the 
Fourier transform of e Ji ^ > ~ J ^ > : 



p(E) = j _^_ e J(t)-J(0)+iEt/h 



(5) 



The function P(E) gives the probability that the envi- 
ronment absorbs a quantity of energy equal to E during 



the tunneling process. 
P{E) reads^i 



Finally the current in terms of 



I 



1 



cRt J _ 



EdEP(eV - E) 



(1 



-0eV\ 



1 



-/3E 



(6) 



where Rt is the tunneling resistance in the absence of 
the environment [i.e. when P(E) = 5(E)]. 

In order to find the behavior of the device we need to 
determine the function P(E) in the presence of the os- 
cillator. Since we will restrict to the quadratic part of 
the Hamiltonian, we can obtain the quantum phase cor- 
relation function J(t) from the classical response func- 
tion of the phase to a current source at the junction. 
Specifically wc add to the total Hamiltonian ([2]) the term 
— Iext(t)4>mt{t)h/ e, where <^i n t(i) is the phase in the inter- 
acting picture. The linear response to the external field 



:t gives 



4>(t) 



dt'x(t-t')I ex t(t')- (7) 
e 



with X (t) = %{[M0),Mt)])6(t). Using 
fluctuation-dissipation theorem one can show that 



J(w) 



-2h 



Im[xH] 
1 - 



the 



(8) 



Since the Hamiltonian is quadratic, x( w ) can be calcu- 
lated from the classical response function (see for instance 
Ref. HI) . The problem is now reduced to the determina- 
tion of the classical retarded response function. This can 
be readily obtained by solving the classical equations of 
motion by Laplace transform. 

We determine now the equations of motion for the 
phase, the position of the oscillator, and the degrees of 
freedom of the environment. To describe the dissipation 
we introduce a bath of harmonic oscillators coupled to the 
mechanical oscillator and a bath of electromagnetic oscil- 
lators coupled to the current through the tunnel junction. 
The coupling of the mechanical and electrical part comes 
from the charge energy Q 2 /2C(x) of the capacitor. The 
equations of motion read: 

mx = -kx -^2k p (x - x p ) + -C'(x)<fi 2 (-) (9) 



C(x)4> 



L n 



+ Ext(t)- 



-l—(<t>e ~ 0, 
J->n 



(10) 

(11) 

(12) 



Here x p and 4> n are the degrees of freedom of the envi- 
ronment, (j> e (t) = eVt/h is due to the constant external 
voltage source and we choose the zero of x as the equi- 
librium position of the oscillator when V = 0. We also 
define C'(x) = dC/dx. 
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For I^t = 0, this system of equations has a stationary 



solution given by x 



C'(x eq )V 2 /2k, = eV/h, 



x. p = x eq and <f> n — 0. This solution is unique for 
a given gauge. We thus expand the equations around 
this solution and consider the quadratic fluctuations of 
the fields only: x(t) = x eq + i{t), x p (t) = x eq + x p (t), 
and <p(t) = eVt/h + <p(t). Let us introduce the Laplace 
transform for the fluctuating fields. For example, for x 
it reads: x(s) = / + °° dte~ st x(t), with Re[s] > and 

x{t) = j+l™^{ds/2m)e st Z(s) with a > 0. The sys- 
tem of differential equations then becomes an algebraic 
system of equations. We are interested in the response 
function. Thus, we can neglect transient terms. Keeping 
only linear terms in the fluctuating fields, and solving 
first for the degrees of freedom of the environment one 
finds: 



[s 2 +j(s)s + ul l ] x 



cvh ~ 



(13) 



.eV 



[Cs + Y(s)] 8 <j> = -C'sx—+I ext (t)-,(U) 

where we have introduced the memory function of the 
oscillator and the impedance of the circuit: 



7(s) 



1 



Z(s/iY 



1 



(Sect. MI A[) . the electromagnetic environment is an LC- 
oscillator. Thus, the system reduces to two coupled os- 
cillators: one mechanical and the other electromagnetic. 
This example is interesting for his simplicity. It clearly 
shows how the mechanical and electric part becomes cou- 
pled and how this coupling appears in the current voltage 
characteristics. By contrast, its experimental realization 
is non-trivial, since it is in general not easy to get rid 
of parasitic capacitances and resistances. An inductance 
and a capacitance are present in the environment for the 
oscillating tunneling junction in the experiment of Ref. 
|33| . The environment forms in this way a radio frequency 
resonator that is used to improve the band-width of the 
detection scheme. However, one should note that ex- 
periment of Ref. [33] is not in the regime of parameters 
considered in the present paper. 

In the second case (Sect. MI Bp . the environment pro- 
vides the dissipation and the conditions for the standard 
Coulomb blockade with the suppression of the current at 
low voltages. This case is more relevant from the experi- 
mental point of view, since high impedance environment 
can be more easily devised. For electronic transport with- 
out mechanical motion an accurate experimental study of 
the effect of a purely resistive environment has been car- 
ried out recently by tuning the environment resistance 
for a given tunnel junction^ 



The two coupled equations (fTBj) and (|14|) describe the re- 
sponse of the phase </> to the external current I ext flowing 
through the tunnel junction. Solving for (f> and letting 
s — > icu + + wc obtain the expression for x(w). Using 
Eq. ([8]), the expression for J(lo) reads: 



47rRe[Z t (a;)] 



1 



1 



(15) 



with 



C' 2 V 2 im/m 



iuiC + Zloj) 



(16) 

The usual result of the Coulomb blockade can be ob- 
tained by setting C = in this expression. Then, only 
the electromagnetic environment would contribute to the 
response function. The presence of the oscillator simply 
modifies the response of the system; for 7 constant, it 
would be equivalent to a RLC circuit in parallel with the 
tunnel junction. 

The current voltage characteristics can now be ob- 
tained by combining Eqs. ©, ©, (fT5|) . and (fTB]) . In the 
following section we consider two representative cases for 
the circuit impedance and discuss the expected conduc- 
tances. 



III. RESULTS IN SPECIFIC CASES 

We consider now two simple cases for the impedance 
of the circuit: Z — itoL and Z = R. In the first case 



A. Inductive case Z{uj) — iLuj 

We consider the case of undamped harmonic oscilla- 
tor. Strictly speaking, the hypothesis of thermal equilib- 
rium for the oscillator and for the electromagnetic modes 
is not valid. We will assume that an infinitesimal dis- 
sipation is present and that the tunneling resistance is 
sufficiently large to let the time for the system to relax 
between two tunneling events. The electromagnetic part 
of the circuit is characterized by the oscillator frequency 
ljlc = 1/VLC. The pure electromagnetic case is dis- 
cussed in the review paper Ref. H^- Due to the coupling 
with the mechanical mode two resonant frequencies ap- 
pear, as can be seen by finding the poles of Z t (u>) [cfr. 

Eq. nr 



(17) 

Here we introduced the "coupling frequency" to 2 = 
C' 2 V 2 /Cm. The real part of Z t reduces to a sum of 
four delta functions 



Rq 



E 

<T=± 



[S(u-u a ) + S(w + u; a )] (18) 



with p± = (E c /h J j±)\w 2 ± - uj 2 m \/\uj 2 ± - uj 2 t \. 

It is convenient to express the coupling of mechanical 
and electronic degrees of freedom in terms of a parameter 
that depends only on the device and docs not depend on 



the voltage bias. The reason is that it is easy to vary 
experimentally the voltage for a given device. A phys- 
ically relevant parameter is the variation of the elastic 
energy E e when the oscillator is displaced by a distance 
Ax = Ec(C /C)/k in response to the Coulomb force 
generated by a single electron on the capacitor. With 
this definition we have 



E e = E 2 C {C /C) 2 /2k and 



lE e {eVf/E%. 

(19) 

For a junction made of a STM tip on top of a suspended 
carbon nanotubc, an electromechanical coupling with the 
radial breathing mode of the tube was reporte d 37 ' 38 . In 
this context, u> m and E c were of the same order of mag- 
nitude (around lOmeV). Thus, we estimate the cou- 
pling parameter E e /E c ~ (xq/L c ) 2 , where xq is the zero- 
point motion of the mode (a fraction of Angstrom) and 
L c = C/C (a fraction of nm). In the following, we as- 
sume a (rather optimistic) ratio of E e /E c < 0.1 

For small value of E e /Ec and uj m 7^ lolc the new 
frequencies u>' m and to' LC are only weakly modified by the 
interaction with respect to their bare values: 



'2/2 1, 1 11 1 2 \ 



J LC 



/u LC = 1 - Wj/(w m - CJ 



LC) 



(20) 
(21) 



(note that uj' m = u + if uj m > lolc and u)' m = oj- if 
u> m < ojlc) and the weights of the poles simplify: 



Pn 



Ec 



< PLC 



Er 



J LC 



LC 



(22) 



The inequality holds if the bare frequencies u> m and ujlc 
are of the same order of magnitude. On the other hand, 
when uj m = lolc 3> ^ii we find cu± = ui m ± u>i/2 and 

P+ = P- = Plc/2- 

^From Eq. (fT5|) . one can obtain J(t) — J(0) = 
EoPJAt) withiM 



ja(t) = COth 



f3ftt) a 



(cos(w CT f) — 1) — i sm(co a t) . (23) 



We consider now the zero temperature limit. In this case, 
the function P(E) obtained from Eq. ((5]) has a simple 
expression in terms of an infinite series: 




FIG. 2: Differential conductance of the device as a function 
of the bias voltage for the case Z = iu>L. The bare oscillating 
frequencies are huj m — 0.5Ec and Hujlc = Ec- The elec- 
tromechanical coupling is E e — (full line) and E e — 0.1 Ec 
(dashed line) (in this case, signatures of mechanical coupling 
are not visible for hui m below 0.5Ec)- 



with a suppression of the conductance at low voltage and 
RtG = 1 for eV S> Ec- The interesting situation is 
when plc is of the order of 1. In this case steps are 
clearly visible for eV - E c (cf Fig. H). If p m ^ 
and u> m < lolc the mechanical oscillator manifests itself 
as a series of additional steps of width fnu' m particularly 
visible at the beginning of each step of the LC circuit. 
In practice the picture is slightly more complex due to 
the dependence on V of the four parameters p a and uj a . 
The renormalization of the frequency and the change of 
p introduces a shift of the steps and a weak smooth In- 
dependence of the conductance plateaux, as shown in 
Fig. [21 The presence of the mechanical oscillator thus 
gives a similar result to what has been predicted for a 
single electron transistoi— with Frank-Condon steps and 
a reduction of the current at low voltage. This result is 
found here in a much simpler structure: a single tunnel 
junction coupled in series with an external inductance. 



P(E) 



^2 Pn(Pm)Pn'(pLc)S(E - nhw' m - n'ftw' LC ) 

n,n'—0 

(24) 

where p n (p) = e 9 p n jn\. 

The differential conductance G = dl/dV can be ob- 
tained from the current ([6]). Neglecting the weak In- 
dependence of P(j and uo a in the derivative we obtain: 



R T G{V) 



E 

hnuj ' m + hn' uj ' L c < e V 



Pn(Pm)Pn'(PLc) ■ (25) 



It leads to a double series of steps spaced by hw m and 
fuoLc- If Pm = 0, one recovers the usual sequence of steps 



B. Ohmic case Z(oj) = R 



Let us now consider the pure Ohmic case for the exter- 
nal impedance Z(u>) = R. The first step is again to find 
the poles of Z t (u). For 7(0;) constant, three poles are 
present: one is pure imaginary, 17^, and the two other 
ones are complex, u>± = ±u>' m + ij' In the relevant 
limit of small coupling u>j <C u> m , jrc an d low mcchani- 
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cal damping 7 < w m , jrc we get: 



IRC 



IRC 



TrC + "1 



1 



'Trc 



7=7 + 



IRC + w m 



(26) 

(27) 
(28) 



The mechanical oscillator acquires a damping even if 
7 = 0. This intrinsic damping is proportional to the cou- 
pling. ^From here, we drop the bare mechanical damping 
since it does not induce new effects compared to those 
already introduced by the presence of R. An important 
parameter for this device is the dimcnsionlcss conduc- 
tance of the tunnel junction g = Rq/R. For g 3> 1 the 
Coulomb blockade is no more observable, since the charge 
can vary continuously. In the absence of the oscillator the 
appearance of the Coulomb blockade is signaled by the 
change of the low-energy behavior of P(E), that goes like 
■^2/g-i 34 Thus for g < 2 the current is suppressed at low 
voltage V < Ec/e. If we consider the interesting case of 
huj m ~ Ec, then the condition j « 1 implies also that 
7RC ^ w ?n • in this limit the real part of Z t is simply the 
sum of three Lorentzian functions: 



MZM\ e c 

5 = -7- 2 1rc) 

+ SL [L ( W _ u > y) + L ( W + u > y } 

7RC 



(29) 



with e = i RC hRC, and L(u,i) = ( 7 /2)/(w 2 + 7 2 /4)/tt. 
Since the coupling is weak, the damping of the me- 
chanical mode is also very small and the Lorentzian 
function can be approximated by a Dirac delta-function 
[L(w,7 — » 0) = S(cj)]. At zero temperature this gives for 
the phase correlation function: 



J(t) - J(0) 



Ec 



jo(lRct) +p m [e~ 



where p m = |^^- and j (r) =2/™ 

in tl 

mechanical oscillator. P(E) thus reads: 



1] (30) 



rdx is the 



correlation function of the circuit in the absence of the 



P (E) = E 



Pn(Pm) 



P [(E-nhiu> m )/e} (31) 



where 



1 f + °° 

P (E) = / ^ L f e iEt/h+{E c /vh'yRc)jo{tlRc) (32) 

27rfr J.^ 

is the function P(E) in the absence of the oscillator. For 
the case at hand of small g the function P (E) is peaked 
at E f=a £c with a width ~ (So /n) \/2g\ \ogg\ that van- 
ishes for g — > 0. Thus, the resulting P(-E) is a sequence 
of peaks shifted by huj m . Figure [3] shows the dependence 




FIG. 3: P{E) is displayed for different values of V. We take 
g = 0.015, hu) m = O.LEc and E e = Q.lEo- 



of P{E) on S and V as obtained by numerical integra- 
tion. For eV/E c < ^/E c /E e (Eqs. JSTJ) and (29])) the 
mechanical coupling is negligible and one recovers the 
usual picture of Coulomb blockade, with a single peak 
at E « Ec- For larger values of eV/Ec new peaks ap- 
pear and their positions slowly drift as a function of V 
towards higher energies. The first peak instead drifts to- 
wards smaller energies, due to the e factor in (|3ip . Since 
Trc reduces for larger values of V, the width of the peaks 
also reduces. 

In Fig. [J] we show the same function, but for a higher 
mechanical frequency Hu> m /Ec = 0.5. It shows the ap- 
pearance of the peaks and their shrinking. 

An important case is g — > 0, for which P (E) = 
S(E — Ec)- This is clearly the most favorable case to ob- 
serve the quantum fluctuations of the oscillator directly 
in the I-V characteristics. Indeed, even in the weak elec- 
tromechanical coupling the sharpness of the peaks gives 
rise to steps in G(V), like in the case of a purely electro- 
magnetic resonator. The conductance has thus the form 
(again neglecting the small voltage dependence of p m and 
huj m in performing the derivative): 



R T G(V) 



E 



Pn{pm) 



(33) 



hnu' m <eV —eEc 



the result is plotted in the insert of Fig. [5] for E e /Ec = 
0.1 and g = 10" 4 . 

For the general case the conductance can be obtained 
numerically by integrating P(E). It is shown for three 
different values of the ratio E e /Ec in Fig. [5] The pres- 
ence of the mechanical coupling introduces oscillations 
of the conductance on the scale eV ~ huj m that are pre- 
cursors of the step behavior in the extreme g = case. 
The small increase of the conductance at low voltage is 
the consequence of the rcnormalisation of Ec due to the 
electromechanical coupling (the e factor). 

A last comment is in order for the resistive case. Even 
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FIG. 4: P(E) is obtained numerically for different values of 
V for g = 0.015, hw m = 0.5E C and E e = 0.1E C . 



\(J = 0.015 No coupling 

"^^^^™ ' " to m = 0.LE c , E e = 0.1£ c 

hu m = Q.5E C , E c = 0.1E C 
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'J — - ha;,,, = 0:2Ec, E c - OAEc 




b = io- 4 || 



2 3 
eV/Ec . 
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FIG. 5: Main plot: G(V) for the same parameters as in 
Fig. @ (dashed line) and Fig. @ (dot-dashed line) and for 
vanishing coupling (full line). Inset: GiV) for g = 10 -4 from 
the analytical expression (|33[1 for hui m = O.liSc, E e = (full 
line) and _B e = 0.1i?c (dashed line). The conductance is much 
more sensitive to the electromechanical coupling in presence 
of a high impedance environment. 



if this system is described by a very similar electric cir- 
cuit to that used for a single electron transistor without 
gate, some differences should not be overlooked. In the 
single electron transistor the two resistance are tunnel 
junctions, while here one is a tunnel junction and one 
a ohmic resistor. A tunnel junction is substantially dif- 
ferent from a ohmic resistor, since in the first case the 
charge can be transferred only by a sudden tunneling 
event, while in the second the charge slowly leaks through 
the resistor. The charge quantization is thus not enforced 
by the Ohmic resistor, and this has consequences also on 
the time dependence of the force acting on the mobile 
part. For these reasons our results complement those 



obtained by others authors for suspended single electron 
transistors^l^Li^s 



C. Range of validity 

Let us return to the comparison of position dependence 
of the tunneling and capacitive terms. Using the explicit 
expression (|28|) for 7', we find that the condition to ne- 
glect the effect of the tunneling force (case -fuc <C w m ) 
is: 



1 RqE c 



lB , c RTE e eW 2 



(34) 



This gives the explicit condition on the voltage in order to 
neglect the position dependent part of the Hamiltonian. 
Formally, it is clear that it can be always fulfilled for Rt 
sufficiently large. 

A second possible source of deviation from the behav- 
ior we found comes from the cubic term of of the inter- 
action energy. We neglected the term (</>) 2 with respect 

to eV4>/h. This is justified when ((</>) 2 ) 1/2 < eV/h. We 
thus obtain the condition of validity of the quadratic the- 
ory at zero temperature: 



— » y/g\]ng\ 



(35) 



IV. JOSEPHSON JUNCTION 

Observing a Joscphson supercurrent through a STM 
tunnel junction has been shown to be feasible.— Par- 
ticularly appealing is the possibility of studying an os- 
cillating Josephson junction. The effect of the position 
dependence of the electronic tunneling amplitude on the 
dynamics of the Josephson junction has been considered 
recently^ However the effect of the position dependence 
of the capacitance (Coulomb interaction) has not been 
explored so far. The theory of Coulomb blockade for 
normal metallic electrodes presented above can be gen- 
eralized with little change to the superconducting case. 
The main difference is that the phase difference is now a 
truly observable quantity and it is conjugated with the 
number of Cooper pairs that cross the junction. 

The Josephson junction is characterized by a critical 
current I c that defines the energy scale Ej = hl c /2e. 
For Ej/Ec <C 1 a tunneling approach can be used to 
obtain the dissipative current of the device. We remind 
the main steps of the procedure^ The current opera- 
tor is Is = I c sin(20) while the Josephson part of the 
Hamiltonian reads Hj = Ej cos(20) (we keep the same 
definition of the phase, the factor two intervenes in the 
usual Josephson relation). 

Calculating the current in the linear response limit we 
get 



Ej 



+00 



at' 



-2i<t>(t') 



(36) 
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FIG. 6: /(V) is plotted for g = 0.065. We distinguish the case 
E e = (full line), E e = 0.1-Ec with hu m = 0.1£c (dashed 
line) and E e — O.lEc with fiw m = 0.5-Ec (dot-dashed line). 
As long as hui m remains smaller than the width of the peak, 
the coupling manifests itself in a broadening of the Coulomb 
peak. For hui m /Ec — 0.5, contributions of the mechanical 
oscillator clearly appear. 



One can thus express the current again in terms of the 
Fourier transform of the phase correlator: 

I(V) = [P s (2eV) - P s (-2eV)] (37) 

where the definition of P S (E) differs from the definition 
of P{E) simply by a factor of 4 in front of J(t): 

/ + OC 
dte 4J(t)-4J(°)+iBt/H (38) 
-oo 

The function J(t) depends on the mechanical response 
function as shown in Section [TTJ The results of Sects. ITT1 
and IIIII can thus be applied to the superconducting case 
to obtain P S (E) and hence I(V). In the Figure [5] we 
show the expected current-voltage characteristics for an 
oscillating Josephson junction in the pure resistive case of 
Sect. IIIII The superconducting device is much more sen- 
sitive to the presence of the mechanical oscillator, since 
the voltage dependence is directly related to P S (E), and 



not to its integral. One can then observe directly the 
Coulomb blockade peaks in the I-V characteristics. 



V. CONCLUSION 

In this work, we considered the current- voltage charac- 
teristics of an oscillating tunnel junction. At low tunnel 
coupling and moderately low voltage bias, we found that 
the coupling between the electronic and mechanic degrees 
of freedom mainly originates from the position depen- 
dence of the junction capacitance. That is, the electrons 
crossing the tunnel junction couple to an electromagnetic 
environment formed by the rest of the circuit which is it- 
self coupled to the mechanical oscillator. In our study, 
we included both the electromagnetic and mechanical en- 
vironments in the Coulomb blockade theory. We found 
that the mechanical properties of the junction result in an 
additional mechanism for the current suppression that is 
similar to the Frank-Condon blockade predicted in an os- 
cillating single-electron transistor. As in the conventional 
Coulomb blockade, similar, though more pronounced, ef- 
fects would also appear in a voltage-biased oscillating 
Josephson junction. 

At larger tunneling rates, two new ingredients should 
be added to our study. First, the electromagnetic and 
mechanical degrees of freedom would not have time to 
relax to equilibrium between two tunneling events. Sec- 
ond, the position dependence of the tunneling matrix el- 
ement may not give anymore negligible effects compared 
to those arising from the position dependence of the ca- 
pacitance. Considering both, it may be interesting to 
investigate whether a shuttling instability takes place in 
a single oscillating tunnel junction. 
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